Exceptional points in the spectra of atoms in external fields 
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We investigate exceptional points, which are branch point singularities of two resonance eigen- 
states, in spectra of the hydrogen atom in crossed external electric and magnetic fields. A procedure 
to systematically search for exceptional points is presented, and their existence is proven. The prop- 
erties of the branch point singularities are discussed with effective low-dimensional matrix models, 
their relation with avoided level crossings is analyzed, and their influence on dipole matrix elements 
and the photoionization cross section is investigated. Furthermore, the rare case of a connection 
between three resonances almost forming a triple-degeneracy in the form of a cubic root branch 
point is discussed. 
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I. INTRODUCTION 
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resonance spectra, namely "exceptional points" [l|,J; 
recently attracted growing attention theoretical^ 



A special type of degeneracy in parameter-dependent 

ional points" [l|,j2|jias 
)n theoretically [J, [ 1 , [ a , 
i,0,[i,Q as well as"experimentally [11 [ul E [ll Qllf ■ 

Exceptional points can be found in systems which depend 
on at least two real valued parameters. They are posi- 
tions in the parameter space at which usually two energy 
eigenvalues pass through a branch point singularity, i.e., 
the two eigenvalues can mathematically be described by 
two branches of the same analytic function and the excep- 
tional point represents the branch point singularity. In 
this case, also the eigenvectors, or wave functions, pass 
through a branch point singularity and, hence, are iden- 
tical. The occurrence of a geometric phase [la] is one of 
the important consequences of exceptional points [ll| . 

Within the framework of the linear Schrodinger equa- 
tion, the existence of resonances, i.e., decaying unbound 
states is important for the appearance of exceptional 
points because the coalescence of two discrete eigenstates 
with identical eigenvectors is not possible in the spectra 
of Hermitian Hamiltonians with potentials, which de- 
scribe bound states. Here, always a set of orthogonal 
eigenstates, which never can become identical, exists. For 
resonances, the situation is different. Their eigenstates 
can become identical. One possibility of obtaiiiing the 
resonances is the complex rotation method [13, llSl. [l9j. 
which leads to non-Hermitian Hamiltonians in the spec- 
tra of which resonances are uncovered as complex eigen- 
values. 

Physical systems in which exceptional points can ap- 
pear have to depend on at least a two-dimensional pa- 
rameter space. If exceptional points are to be observ- 
able it must be possible to adjust these parameters in 
a sufficiently wide range of the parameter space. Fur- 
thermore, the complex energy eigenvalues, typically the 
positions (frequencies or energies) and widths of reso- 
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nances, must be accessible with a high precision. Exam- 
ples are discussed, e. g., for complex atoms in laser fields 
[2(j | , a double S well 21 1 , the scattering of a beam of par- 
ticles by a double barrier potential [22l |. non-Hermitian 
Bose-Hubbard models [23|, or models used in nuclear 
physics [2^ . The resonant behavior of atom waves in 
optical lattices (25|l shows structures originating from ex- 
ceptional points [5], and they can be found in nonlin- 
ear quantum systems. The stationary solutions of the 
Gross-Pitaevskii equation describing Bose-Einstein con- 
densates exhibits a coalescence of two states due to the 
nonlinearity of the equation, which turns out to be a 
branch point singularity of the energy eigenvalues and 
wave functions [26l[27j. However, the phenomenon of ex- 
ceptional points in physics is not restricted to quantum 
mechanics. Acoustic modes in absorptive media 28] rep- 
resent a mechanical system, in which branch point singu- 
larities appear. Manifestations of exceptional points can 
also be seen in optical devices [H [H, jsij, [12] . The most 
detailed experimental analysis of exceptional points has 
been carried out for the resonances of microwave cavities 
[Tol IT^ . IT3 ] , which open the possibility of studying the 
properties of the complex resonance frequencies and the 
wave functions. In particular, the geometric phase has 
been demonstrated experimentally [111 . [l3[ . 

Atoms in static external electric and magnetic fields 
are fundamental physical systems. As real quantum 
systems they are accessible both to experimental and 
theoretical methods and have, e.g., very successfully 
been used for comparisons with semiclassical theories 
[3^ . m, [H, HI]. They provide a rich variety of physi- 
cal phenomena such as Ericson fluctuations, which were 
discovered both in numerical studies [stI , [ssj and experi- 
ments [39] . Recently, the existence of exceptional points 
was reported and a method to detect them in an exper- 
iment with atoms was proposed [40| . Atoms in external 
fields are well suited to study the signatures originating 
from exceptional points in quantum spectra. The electric 
and magnetic field strengths span the two-dimensional 
parameter space required to provide the branch points. 

It is the purpose of this paper to show how excep- 
tional points can be detected systematically in spectra 
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of the hydrogen atom in external fields and to discuss 
their properties in detail. Numerically exact resonance 
energies and eigenstates are determined by the complex 
rotation method (13, HI, The permutation behav- 
ior of the resonances for closed parameter space loops is 
used to unambiguously find and verify the branch point 
singularities. Examples for exceptional points found in 
the system are presented. The shapes of the loops of 
the complex eigenvalues representing the resonances are 
explained with a two-dimensional matrix model. It is 
shown that single dipole matrix elements of two states 
coalescing at an exceptional point diverge whereas this 
behavior does not carry over to the observable photoion- 
ization cross section. 

Besides the typical case, where an exceptional point 
consists of two resonances forming a square root branch 
point, higher degeneracies are possible. In this paper 
the rare case of a structure with three resonances almost 
forming a triple-coalescence in the form of a cubic root 
branch point is presented. It is shown how exceptional 
points are related to avoided crossings of the energies or 
widths of the resonances. 

The paper is organized as follows. Sec. [TTl summarizes 
the most important characteristics of exceptional points 
and demonstrates them with the help of an illustrative 
example. The Hamiltonian of the crossed-fields hydrogen 
atom is introduced in Sec. Illli and the method to obtain 
the resonance eigenstates is presented. A procedure to 
find exceptional points in spectra of atoms in external 
fields is presented and applied in Sec. IIVI in which also 
examples of branch point singularities are given. A two- 
dimensional matrix model used to explain the paths of 
the eigenvalues for parameter space loops is introduced 
and applied in Sec.|Vl The influence of exceptional points 
on dipole matrix elements and the photoionization cross 
section is discussed in Sec. |Vll Three resonances al- 
most forming a triple-coalescence are presented in Sec. 
IVIIl and the connection between exceptional points and 
avoided level crossings is investigated in Sec. IVIIII Con- 
clusions are drawn in Sec. IIXI 



II. EXCEPTIONAL POINTS 
A. Important properties 

The typical case of an exceptional point is a position in 
a two-dimensional parameter space at which two eigen- 
values pass through a branch point singularity. This 
behavior has important consequences for the associated 
eigenvalues and eigenvectors. 

If one encircles the exceptional point in the parameter 
space, a typical permutation behavior of the eigenvalues 
can be observed [3] ■ The two eigenvalues which represent 
the two branches of one analytic function with the singu- 
larity are interchanged after one circle around the excep- 
tional point, whereas all further eigenvalues do not un- 
dergo a permutation. Away from the exceptional points 



the branching eigenvalues are different and each of them 
belongs to a distinct eigenvector. Accordingly, these 
eigenvectors undergo the same permutation as the eigen- 
values, and at the exceptional points they likewise pass 
through a branch point singularity 0, i.e., there is only 
one linearly independent eigenvector for the two degener- 
ate eigenvalues. In a matrix representation, the eigenval- 
ues and eigenvectors form a normal block [4l|, \^ , which 
distinguishes them from ordinary degeneracies where two 
eigenvalues belonging to two different analytic functions 
have the same value. Alongside the permutation of the 
eigenvectors a geometric phase appears. In the most com- 
mon physical situation of a square root branch point re- 
sulting from a complex symmetric matrix, the geometric 
phase becomes manifest in the change in sign of one of 
the two eigenvectors ^)] and can be written in the form 

r ^ circle r ^ 

[Xi,X2\ ~> [X2, ~Xi\ , 

where Xi and X2 are the eigenvectors permuted during 
the loop around the exceptional point. 

It should be mentioned that the occurrence of the phe- 
nomenon is not restricted to two dimensions. In higher- 
dimensional parameter spaces of complex matrices, an 
exceptional "point" is always an object of co-dimension 
2 [42]. That is, in two dimensions an exceptional point 
indeed appears as a point, whereas in three dimensions 
it is a one-dimensional object, or line, and in a four- 
dimensional parameter space it has the form of a two- 
dimensional "exceptional surface." 

B. A simple example in a non-Hermitian linear 
map 

One of the simplest examples in which an exceptional 
point occurs, and which helps to understand many effects 
related to it, is described by the two-dimensional matrix 

a 

M(«)^(^_^) (1) 

with the complex parameter k. The two eigenvalues of 
the matrix are given by 

Ai = v/TT^, A2 = -Vl + n^ (2) 

and are obviously two branches of the same analytic func- 
tion in K. There are two exceptional points in the system, 
which appear at the complex conjugate values k± = ±i 
as can easily be seen. The eigenvectors which belong to 
the two eigenvalues are 

They also depend on the parameter k and pass through a 
branch point singularity at the exceptional points k± = 
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(a) 1.00006 



1.00000 



0.99994 



(b) 



-0.00006 



0.00006 



Re(K) 




-0.012 



-0.012 



0.012 



of view of exceptional points. The electric field is as- 
sumed to point in the x-direction and the magnetic field 
is orientated along the z-axis. The Hamiltonian without 
relativistic corrections and finite nuclear mass effects [i^l 
reads in atomic Hartree units 



H 



\p'-\ + l^L, + ^7'(a^' + y^) + , (6) 



FIG. 1: (Color online) (a) Circle in the parameter space k 
with the exceptional point k+ = i as center point for the 
simple model ([T}. (b) Eigenvalues Ai,2 calculated for the pa- 
rameter values from (a) indicated by red squares and green 
diamonds, respectively. In this case each of the two eigenval- 
ues traverses a semicircle. In (a) and (b) the filled symbols 
represent the first parameter value kq and the corresponding 
eigenvalues Ai,2(ko), respectively. The arrows point in the 
direction of progression. 



±i, where the only linearly independent eigenvector reads 
a;(±i) = (Ti,l). ... 

The branch point singularity leads to a characteristic 
behavior of the corresponding eigenvalues under changes 
of the parameters. If one chooses a closed loop in the 
parameter space and calculates the eigenvalues for a set 
of parameters on this loop, the permutation of the eigen- 
values can be seen by plotting their paths in the complex 
energy plane. The situation is illustrated in Fig. [TJ Here, 
a circle K{ip) = i + ge^"^ around the singularity k+ = i is 
traversed in the parameter space, which is shown in Fig. 
[TJa). After one revolution, the first eigenvalue marked 
by red squares has traveled to the starting point of the 
second marked by green diamonds, and vice versa. As a 
consequence, the path of each eigenvalue is not closed if 
one traversal of the loop in the parameter space is per- 
formed, but the path is closed if the parameter space loop 
is traversed twice. Of course, for the simple model it is 
also possible to demonstrate the half-circle structure an- 
alytically. If the parameter space curve described above 
is applied to the eigenvalues ^ one obtains, for <C 2, 
the expansion 



Ai,2 = ±Vl + (i + eci^)2 = ±Vec''^^ V2i + pei 



where p is the kinetic momentum of the electron, r its 
distance from the nucleus, and Lz the z-component of 
the angular momentum. The Hamiltonian contains the 
Coulomb potential oc 1/r, the paramagnetic term oc 7^2, 
the diamagnetic term oc 7^(2;^ + y"^), and the potential 
due to the external electric field oc fx. 

For the numerical calculation of the energy eigenvalues 
in a matrix representation of the Hamiltonian ([6]) it is 
advantageous to transform the Schrodinger equation into 
dilated semiparabolic coordinates. 
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U = -\Jr + z , V = -\Jr — z 
b 



If = arctan — (7) 

X 



with r = y'' x"^ + y'^ + z^. For calculating the resonances 
of the system the complex rotation method [13, [H, [3 
is used as a powerful tool |33, HI, 113 • The complex 
scaling is included via the dilation parameter &, where 
the replacement 

= \b^\c'^ (8) 
entails the complex rotation of the position vector 



e r 



(9) 



±v/2^)e'"/*e 



i7r/4 i^/2 



(4) 



and leads to the complex scaled and regularized 
Schrodinger equation 



- 2Ho + Ibl^'^S {^^' + i^') i^ - i|6|«e'4.^2 

X fi'^iy^ + v^) - 2\b\^(^^'^ f^lv (/i^ v"^) cos 

+ 4|6| V« + {^? + v^) je-'^V = 2\b\''E (/i^ + ^ 

(10a) 



or 

Ai = v^e'(^/4+^/2) , A2 = v/2^c'(5V4+^/2) ^ 

which reproduces the half-circles shown in Fig. [T] for a 
full parameter space loop ip = Q . . .27: correctly. 

III. HAMILTONIAN AND MATRIX 
REPRESENTATION 

In this paper the hydrogen atom in static crossed elec- 
tric and magnetic fields is investigated from the point 



with the term 

i^o = -^(A^ + A,) + ^ (m' + , (10b) 

where 

1 d d I 
^'^l>d-g'd-g^7W^ ^^elM,.}. (10c) 

Due to the harmonic oscillator structure of Hq a well 
suited complete basis set is given by the states 

\n^,nu,m) = \nf^,m) ®\ny,m) , (11) 
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where each of |np,m) and |ny,m) represents an eigen- 
state of the two commuting operators 

TV ~ a^ai + aja2 , (12a) 
L ^ i{aia^ - a2) ^ iqiP2 - q2Pi) (12b) 

of the two-dimensional isotropic harmonic oscillator with 
common eigenvalue m of L. The operators Ui and aj' 
are the familiar ladder operators of the one-dimensional 
harmonic oscillator. 

The matrix representation of the Schrodinger equation 
pOap is non-Hermitian and has the form 

A{^,m^2b^EC^ , (13) 

where A(7, /) is a complex symmetric matrix and C is 
real symmetric positive definite. Above the ionization 
threshold, resonances are uncovered as complex energy 
eigenvalues E, where the real and imaginary parts repre- 
sent the positions and the widths T = — 2Im(i?), respec- 
tively. 

The Hamiltonian has two constants of motion, namely 
the energy and the parity with respect to the {z = 0)- 
plane. The latter symmetry opens the possibility of clas- 
sifying eigenstates by their z-parity and to consider the 
associated subspaces separately. The examples discussed 
in this paper are given for even z-parity. 

The computation of the eigenvalues was performed us- 
ing the ARPACK library [4^. For typical calculations, 
the number of basis states was on the order of 10,000 
to 12,300. The exact determination of the positions of 
exceptional points up to four valid digits often required 
larger matrices with up to 17,300 states. If the matrices 
A and C are built up appropriately, they possess a band 
structure, which was exploited in the numerical diago- 
nalizations. 



IV. EXCEPTIONAL POINTS IN SPECTRA OF 
THE HYDROGEN ATOM IN EXTERNAL FIELDS 

A. Identification of exceptional points 

The Hamiltonian in the Schrodinger equation IjlOap is 
non-Hermitian. Two parameters, viz. the strengths of 
the external electric and magnetic fields, are available 
to influence the positions of the resonances and thus it 
should be possible to produce degeneracies of the com- 
plex resonance energies. Exceptional points do exist in 
atomic spectra if the fields can be chosen in such a way 
that a coalescence of two states occurs. The crossed-fields 
hydrogen system fulfills all necessary conditions for the 
appearance of exceptional points, however, one has to 
find them in the spectrum to really prove their existence. 

A successful procedure to search exceptional points can 
be found if one exploits their properties. The permuta- 
tion of the two eigenvalues involved in the singularity 
provides a clear signature which can be used to detect 
exceptional points. A good choice for a closed loop is a 



"circle" in the parameter space of the two field strengths 
with a radius S < 1 chosen relative to a center (70, /o), 

7(v) = 7o(l + (^cosv?) , f{(p) ^ foil + S sin (p) . (14) 

It opens the possibility of scanning a larger area of the 
parameter space, namely the complete circular area, at 
once. A fundamental advantage of the method is that 
it allows for automatizing the procedure up to a certain 
extent. If the steps on the circle are chosen small enough, 
the resonances of two consecutive steps can be assigned to 
each other unambiguously, and their motion in the com- 
plex energy plane can be traced. Eigenvalues which do 
not return to their starting point once the circle in param- 
eter space is closed but are interchanged with a further 
resonance are a proof for the existence of an exceptional 
point. The exact position of the degeneracy can be deter- 
mined by minimizing the distance of the two eigenvalues. 
After the degeneracy has been found numerically, a last 
circle with a small radius (typically S « 10~^^) around 
the parameter point at which the degeneracy occurs is 
used to decide whether or not the branch point singular- 
ity structure is present and the degeneracy found is an 
exceptional point. 

The geometric phase appearing with exceptional points 
is accessible through the eigenvectors representing the 
resonances in the numerical calculations, and provides a 
further possibility of v erify ing their existence [40]. As 
was also shown in Ref. [401] the permutation of the com- 
plex resonance energies opens the possibility of detecting 
exceptional points in experiments with atoms once the 
complex energies have been extracted from the photoion- 
ization cross section. 



B. Examples 

With the method described above, exceptional points 
have been found in spectra of the hydrogen atom in static 
external fields Table [J lists 17 examples. In the cal- 
culations the relative difference — i?2 | of the two 
eigenvalues could be reduced down to 10^ However, 
this is only the result for a single matrix representation. 
What is more crucial is the influence of the complex ro- 
tation on the matrix with finite size. Using up to 17,300 
states the convergence of typically three to four valid 
digits in the parameters as well as in the energies can be 
achieved. The convergence was checked with the stabil- 
ity of the results against changes in the matrix size and 
the complex parameter b. 

As an example Fig.[2]shows a typical result obtained in 
a numerical calculation for the exceptional point labeled 
14 in Table [J The squares and the diamonds represent 
each of the eigenvalues at different field strengths. In this 
example, using 20 steps on the circle in the parameter 
space has been sufficient to obtain a clear signature of 
the branch point singularity. The "radius" of the circle 
according to Eq. ^1)1 was 6 = 0.01. 
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TABLE I: Examples for exceptional points in spectra of the 
hydrogen atom in crossed magnetic (7) and electric (/) fields. 
All values in atomic units. The numbers are used as labels to 
identify the exceptional points. 







f 

J 




Im(£') 


1 


002335 


0001177 


—0.01767 


—0 000103 


2 


0.002575 


0.000117114 


—0.015067 


—0.0000823 


3 


002752 


0001298 


—0 015714 


—0 00022637 


4 


0030152 


0001231 


—0 01209 


—0 000099 

yjtyj\jyj\j\j\j 


5 


003045 


0001332 


—0 015812 


—0 0001896 

yjtyj\jyjA-^\jyj 


6 


0.0030460 


0.000127302 


—0.017624 


—0.000087 


7 


0037Q1 5 




—0 fll 24fl 


—0 nnni 64 

yj •\jyJ\j -L U4: 


8 


0.004604 


0.0002177 


-0.022135 


-0.00006878 


9 


0.004714 


0.00021529 


-0.01394 


-0.00010 


10 


0.00483 


0.000213 


-0.01255 


-0.00030 


11 


0.00529 


0.0002011 


-0.0150 


-0.000136 


12 


0.00537 


0.000214 


-0.01884 


-0.0000679 


13 


0.005388 


0.0002619 


-0.02360 


-0.00015 


14 


0.00572 


0.000256 


-0.01984 


-0.000258 


15 


0.00611 


0.000256 


-0.01593 


-0.00024 


16 


0.00615 


0.000265 


-0.0158 


-0.000374 


17 


0.00776 


0.000301 


-0.0179 


-0.000756 



-0.00010 
-0.00015 
-0.00020 
-0.00025 
-0.00030 
-0.00035 



-0.00040 

-0.0200 



0.00026 



0.00025 



□ J 



□ □ 



0.0056 0.0058 
y [a. u.] 



o 



o 
o 



o o o 



-0.0199 -0.0198 

Re(£) [a. u.] 



-0.0197 



-0.0196 



FIG. 2: (Color online) Paths of the two eigenvalues (repre- 
sented by squares and diamonds, respectively) which degen- 
erate at the exceptional point labeled 14 in Table |l] in the 
complex energy plane [401]. Each point of one eigenvalue be- 
longs to a different set of parameters. The path in the field 
strength parameter space is a circle defined in Eq. (|14p with 
S = 0.01 (see inset). The initial set of parameters and the 
corresponding eigenvalues are represented by filled starting 
points. The arrows indicate the direction of progression. The 
filled triangle marks the position of the exceptional point in 
the parameter space and the corresponding complex energy 
of the degenerate resonances. 



V. DESCRIPTION OF THE RESONANCE 
ENERGIES IN THE VICINITY OF 
EXCEPTIONAL POINTS 

A. Effective two-dimensional matrix 

The most common case of an exceptional point consists 
of two resonances forming a square root branch point, 
whose effects on the spectrum appear in a close vicin- 
ity of the critical parameter values. There, only the two 
resonances involved in the branch point singularity are 
important if they are sufficiently separated from all fur- 
ther resonances, which in general can be realized since 
only a narrow region of the complex energy plane around 
two almost degenerate eigenvalues must be taken into ac- 
count. In this case, one can restrict the discussion to the 
subspace spanned by the two relevant eigenvectors close 
to the exceptional point, which leads to effective two- 
dimensional matrix models. Indeed, two-dimensional 
models with only one complex parameter similar to the 
model introduced in Sec. Ill Bl yield a good description of 
the two complex eigenvalues in the vicinity of exceptional 
points [l| , and they can also provide good results for the 
resonances of the hydrogen atom in external fields. How- 
ever, for some effects the actual matrix structure, which 
includes the two real field strengths 7 and /, has to be 
taken into account. We introduce a model adequate to 
describe the physical crossed-fields system. It is given by 
a matrix, whose elements have the form 



,(0) 



+ «.?(7 - 70) + c^i'if - /o) , {1, 2} . 

(15) 

This two-dimensional matrix is still a very simple model 
because it includes only the two states merging at the 
exceptional point and ignores couplings to other levels. 
Furthermore, a linear dependence of the matrix elements 
on the two field strengths is assumed, which is certainly 
only true for small distances to the center point (70, /o)- 
In contrast with the simple model used in Sec. Ill Bl it in- 
cludes two real parameters with complex prefactors a\^"^ , 

a\j \ a\P and correctly reproduces the matrix shape of 
the full Hamiltonian to lowest order. In a power series 
expansion its eigenvalues fulfill the relations 

Ai + A2 = tr(M) = A + 5(7 - 70) + C{f - fo) , 

(16a) 

(Ai - Aa)^ = tr(M)2 - 4det(M) = D + E{-f ~ 70) 
+ F{f - fo) + G(7 - 70)' 
+ ^(7-7o)(/-/o)+/(/-/o)' (16b) 

with new coefficients A, B, C, D, E, F, G, H, and /. 
Since the eigenvalues do not change under a similarity 
transformation of the matrix M and the explicit choice 
of the matrix is not relevant, the representation (|16ap 
and (|16bp is more suitable than Eq. (jlSp . in which more 
coefficients appear. The coefficients can be determined 
by a fit of the exact quantum energies to Eqs. p6ap and 



,(/)/ 
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(a) -0.00005 




(b) -0.00015 



-0.00050 I 

-0.0203 Re(£) [a.u.l -0.0195 

(c) -0.000245 



¥ 

-0.000270 



-0.00040 I 

-0.0200 Re(£) [a.u.] -0.0197 




-0.01985 Re(£:) [a.u.] -0.01983 

FIG. 3: (Color online) Paths of the two eigenvalues consid- 
ered in Fig. [2] for different radii of the parameter space circle 
(|14[) . where the center (70, /o) = (7'^^\/'^^') was always 
chosen to be exactly the exceptional point, (a) 5 = 3 x 10""^, 
(b) 5 = 10"^ (c) 5 = 10"". The position of the reso- 
nances at (7'^^\/'^^') is marked by a triangle in each fig- 
ure. The squares and diamonds represent the exact quantum 
resonances. The lines represent the eigenvalues of the two- 
dimensional matrix model psp . whose coefBcients have been 
fitted to the numerical results of the exact quantum calcula- 
tions. 



(jl6bp in a region of the parameter space in which only two 
resonances are relevant. A fit for 6 different parameter 
sets yields the 9 coefficients yl to / and, thus, determines 
completely the two eigenvalues of the model p3|) in de- 
pendence of the deviations from the center point (70, fo). 
Six differences Ai — A2 are required to determine D, E, 
F, G, H, and /. Three of the same parameter sets can 
be used to determine A, B, and C from the sum Ai -|- A2. 



B. Shapes of the eigenvalue loops 

As can be seen in Fig. [5] the shape of the paths covered 
by the two energy eigenvalues differs considerably from 
the semicircle form obtained for a small loop in the sim- 
ple two-dimensional matrix model in Fig. [ijb). A good 
description of the complicated behavior is possible with 
the model (|16ap . (jl6bp . The results of a more detailed 
investigation of the phenomenon with the matrix model 
are shown in Fig. [3] A circle around the exceptional 
point, which is always located exactly at the center of 
each figure, is performed for three different radii. For 
radii 5 = 0.03 and 6 = 0.01 (see Figs. Efa) and (b)) 
the complicated structure already known from Fig. [2ta) 
appears. The deformations of the eigenvalue paths can 
be reproduced with the matrix p5p . The lines in Figs. 
[3]Ja) and (b) represent the positions of the two model 
eigenvalues Ai and A2 for the same parameter space cir- 
cle which was used for the exact quantum resonances. 
The very good agreement demonstrates that it is possi- 
ble to describe the local structure of the resonances at an 



exceptional point with a simple two-dimensional model 
that ignores the influence of further resonances even if 
complicated structures in the eigenvalue paths appear. 

Fig. [3fc) shows a circle around the same exceptional 
point for the much smaller radius 5 = 10~^, where the 
shape of the paths becomes more similar to the semicircle 
known from Fig. [T]^b) . The exceptional point marked by 
the triangle in Fig. [Jfc) now is located at the center of 
the enclosing eigenvalue trajectories. However, the loop 
still is not a perfect semicircle. This effect is a result of 
the dependence on two real parameters with two com- 
plex prefactors as introduced in the model While 
the models using one complex parameter (cf. Sec. IIIB[) 
lead to a perfect semicircle, as was demonstrated with 
the power series expansion ([5]), this is not the case for 
the description (fT5|) . A short calculation shows that a 
fractional power series expansion similar to Eq. ([4]) for 
the model (|16ap . (|16b[) to lowest order reads 




if |y| < \U\ 

\i\V\ > \u\ 
(17) 

with U = E-ia+ Ffo, and V = E-fo - Ffo- The second 
term under the second square root is large enough to have 
a considerable influence and leads to the modulation of 
the "radius" during the traversal of the semicircle evident 
in Fig. ^c). Again, the model fitted to the numerical 
data, and rendered by the continuous lines in the figure, 
perfectly reproduces the exact behavior. 

Intersections of an eigenvalue path with itself are ob- 
served, which means that the eigenvalue can have the 
same complex energy for two different parameter sets. 
In Fig. ^a) one can even see that the position of the ex- 
ceptional point lies outside the area enclosed by the two 
eigenvalue paths. This is possible if one of the two eigen- 
values (in this example obviously the resonance denoted 
by red squares) has the same complex energy as the one 
at the exceptional point for a second parameter set. The 
crossing of an eigenvalue path with the position of an ex- 
ceptional point in the complex energy plane is shown in 
Fig. m The parameter space circle labeled 1 in Fig. [H^a) 
is chosen such that it passes directly through the point 
(71, /i) at which one of the resonances returns to its ini- 
tial position at the exceptional point as can be seen in 
Fig.m^b). Additionally, the straight line labeled 2 in Fig. 
|4ja) is traversed. For this line one observes that both 
resonances leave the position of the exceptional point. 
The resonance denoted by the solid red line returns to 
its original position on the same path as can be seen in 
the magnified illustration in Fig.|4l^c). 



VI. DIPOLE MATRIX ELEMENTS AND THE 
PHOTOIONIZATION CROSS SECTION AT 
EXCEPTIONAL POINTS 

In the correct inner product, which has to be used for 
the complex rotated states, no complex conjugation of 
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(a) 0.000259 




(b) -0.0001 
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0.000251 

0.00563 Y [a.u.J 0.00582 

(c) -0.000257 



-0.0004 



and the photoionization cross section 




1 



-0.0200 Re(£) [a.u.] -0.0196 



-0.000263 

-0.019861 



«+ 



Re(£) [a.u.] 



-0.019835 



FIG. 4: (Color online) (a) A circle in the parameter space (la- 
beled 1) is chosen such that it passes through the parameter 
set (71 , /i ) at which one of the energy eigenvalues is identi- 
cal with its position at the exceptional point. A straight line 
(labeled 2) connects the parameters of the exceptional point 
and the parameter set (71, /i). (b) Paths of the two numer- 
ically exact complex energies for the circle and the straight 
line described in (a), (c) Magnification of the path of the 
eigenvalue denoted by the red solid line in (b) for the param- 
eter line labeled 2 in (a). The eigenvalue departs from the 
point of degeneracy (circles) and returns to its initial value 
(plus symbols). 



the factors e'^ (cf. Eq. ([9])) may be performed 17]. Only 
the intrinsic complex parts of the wave functions are com- 
plex conjugated. An important consequence is that reso- 
nance wave functions normalized with respect to this in- 
ner product diverge at the branch points, which is typical 
of exceptional points [l[ . In the simple two-dimensional 
matrix model introduced in Sec. Ill Bl that behavior can 
be observed directly for the normalized eigenvectors 



Xi,2(k) 



^i^VTT^) 



2 I It vT 



(18) 

at the exceptional points k = ±i. 

One may wonder whether or not the diverging behavior 
carries over to measurable physical quantities. In partic- 
ular, the photoionization cross section is important for 
the observation of resonances in experiments. For ex- 
ample, in Ref. [i^l it was proposed to measure the pho- 
toionization cross section for parameter sets located on a 
closed curve and to extract the complex energies of the 
resonances from the cross section with the harmonic in- 
version method. This procedure allows for searching the 
permutation behavior in experimental data. 

As presented by Rescigno and McKoy [4^] the dipole 
matrix elements 




can straightforwardly be calculated once the energies ob- 
tained from the complex rotated Hamiltonian and the 
corresponding eigenvectors are at hand. In Eq. (|20|) D is 
the dipole operator in atomic units for a given direction of 
polarization and aoim is the fine-structure constant. The 
bound state, whose energy Eq is supposed to be known, 
is represented by ^o- The ionized states are labeled ^j, 
where the rotated eigenfunctions calculated by the com- 
plex rotation method are used. In converged spectra, the 
result is independent of 9. 

Due to the dependence of Eqs. (fT9| and (|20p on the 



resonance wave functions ^'j a remarkable behavior at 
exceptional points occurs. The diverging behavior of the 
wave functions must carry over to the dipole matrix el- 
ements and, indeed, the numerical results show that the 
resonance wave functions obtained by the complex ro- 
tation method lead to diverging dipole matrix elements. 
This is, however, not an observable physical property be- 
cause the single dipole matrix elements of the two iden- 
tical wave functions at an exceptional point are not ac- 
cessible. In particular, the photoionization cross section 
behaves regularly and does not diverge at an exceptional 
point as will be shown below. 

Again, a two-dimensional matrix model helps in un- 
derstanding these effects. To keep the discussion trans- 
parent, the simplest example, namely the symmetric ma- 
trix model introduced in Sec. Ill Bl is used, the results, 
however, are valid for all complex symmetric matrices. 
Note that in particular there is no difference visible in 
the behavior of the eigenvalues between the lowest-order 
power series expansions (|1]) and ([T7)) if only the distance 
d (or g) is varied and the angle ip is kept constant. The 
normalized eigenvectors, which correspond to the reso- 
nance wave functions '^j, in this model are given by Eq. 
P^ . The form of the dipole matrix elements is 
represented by a product Pi, 2 = y ■ Xi_2 of the eigen- 
vectors Xi^2 with an arbitrary vector y — (?/i,?/2)- If the 
squares of the dipole matrix elements, which are required 
for the photoionization cross section, are calculated for 
a small complex deviation S from one of the two excep- 
tional points, K ~ i + 6, a fractional power series expan- 
sion shows that the single contributions 



,i37r/4 



+ 77(^1+2/2) 



2V2 2 



,i57r/4 



±^{yl-(!>my2^yl)^5 + 0{5) (21) 

diverge as 1/ \/6. It is interesting to note that the sum of 
both contributions always has the exact value 



P,{E) = {■^o\D{r)\^,{e)) 



(19) 



{y-xif + {y-X2f^yl+yl , 



(22) 
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independent of the parameter k, i.e., of the presence of 
the exceptional point, and of the matrix used. What is 
more interesting, however, is the sum 



Ai 



E 



(23) 



which describes the contribution of the two resonances 
to the photoionization cross section with the eigenvalues 
and a real variable E representing the energy. Here, 
one can also look at the contributions of the single eigen- 
values, and obtains in a fractional power series expansion 
around the branch point 



(y • 3^1,2)^ 



i77r/4 



±- 



fi{E,y) 



K2-E 2V2 E^r^ 

±f2{E,y)Vs + h{E,y)6 + OiS''/^) (24) 

with rather complicated expressions fi{E,y), f2{E,y), 
and f3{E,y) which do not depend on S. These parts 
diverge, however, ai and a2 alone are not observable. 
The sum contributes to the photoionization cross section, 
and, in particular, at the exceptional point [S = 0) the 
two resonances overlap. For the sum one finds 



a = 2h{E,y) + 2hiE,y)S + OiS^) 



(25) 



that is, the photoionization cross section converges lin- 
early to a constant value at the branch point. 

Numerical calculations for the hydrogen spectra dis- 
cussed in this paper demonstrate the applicability of the 
simple model. For this purpose, the square modulus of 
the dipole matrix elements (fTO)) and the photoionization 
cross section (pn|) are calculated on a straight line of the 
form (jl4p for a constant angle (p and variable distance 
6 from the exceptional point (70, /o) = (7^^^\ /^^^^)- 
The results have been verified for different angles ip. Fig. 
[5] shows the squares of the two isolated dipole matrix 
elements for the resonances which form the exceptional 
point labeled 12 in Table HI Both matrix elements be- 
have as predicted by the approximation (j2ip of the sim- 
ple two-dimensional model. The single terms and 
diverge in the form of a reciprocal square root of 
5 which is shown by a fit of the numerical results to the 
function 



PatiS) 



+ b 



(26) 



An excellent agreement can be observed. Additionally, 
the real and imaginary parts of the sum + P| are plot- 
ted. As expected from Eq. ((22|) for the two-dimensional 
model, this sum has a constant value. 

The photoionization cross section (7{Re{E^^^'>)) evalu- 
ated at the real part of the energy at the excep- 
tional point according to Eq. (pO)) is shown in Fig. [6] as a 
function of the distance parameter 6. As can be seen di- 
rectly in the figure the numerical data points (red points) 
converge linearly to a constant value for S 0. The red 
line represents a fit to the function 



ffit 



(Re(i?('^P))) = a- 
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FIG. 5: (Color online) Squares P of the dipole matrix el- 
ements between a wave function corresponding to a bound 
state and the two states evolving from the exceptional point 
labeled 12 in Table |I] at field strengths according to the 
form (I14p with the constant angle (p = 0.7 and (70, /o) = 
(7'^^-', Z'^^-'), S is the distance from the critical values of the 
exceptional point. For both states the square modulus di- 
verges in the form of a reciprocal square root of the distance. 
A fit of the data points to a function of the form a/^/S + b, 
which is expected from Eq. (1211) . is shown by the solid red 
and dashed green lines. The real and the imaginary part of 
the sum Pi + P| is also drawn. 




0.000 0.001 0.002 0.003 0.004 0.005 

8 

FIG. 6; (Color online) Photoionization cross section 
a{Re{E'^'^'^'>)) with the bound state energy Eo = -0.125 
(n = 2, p-orbital) as a function of the distance S from the 
critical field values of the exceptional point. The numerical 
data (dots) converge linearly to a constant value as is ex- 
pected from Eq. (|25|) . Again, the line (|14|l with constant 
angle ip = 0.7 and (70, /o) = {'y'-^^K f'-^^'') is used. 



(27) 



whose form is expected from the power series expan- 
sion (|25p of the two-dimensional model. The comparison 
shows an excellent agreement. 
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VII. STRUCTURES WITH THREE 
RESONANCES 



Beyond the typical square root branch point behavior 
studied in most physical examples, higher branch points 
connecting more than two eigenvalues are possible. In 
particular, the possibility of a cubic root branch point 
in a three-dimensional symmetric matrix is the topic of 
current studies [13] ■ For a complex symmetric matrix 
a coalescence of N levels requires {N^ + N — 2)/2 real 
parameters (47j . Thus, for N = 3 five parameters are 
necessary. For the hydrogen atom this means that addi- 
tionally to the two field strengths 7 and / three further 
parameters have to be introduced. Locating exceptional 
points in a five-dimensional parameter space is expected 
to be a laborious task considering the difficulties encoun- 
tered already for two parameters. However, as will be 
shown below, a combination of three resonances strongly 
related to a cubic root branch point is observable in spec- 
tra of the crossed-fields hydrogen atom. In particular 
the permutation behavior of a cubic root branch point, 
where three resonances are permuted and three circles in 
the parameter space are required to restore the original 
situation, can be found. 

An example where three resonances come into play is 
given in Fig. [T] In Fig. [TJa) obviously a permutation 
of three eigenvalues indicated by points with three dif- 
ferent symbols and colors can be observed for a closed 
loop in the parameter space. From that behavior one 
might assume that indeed an exceptional point consist- 
ing of three resonances, which form a cubic root branch 
point singularity, was detected. A more detailed analy- 
sis shows, however, that this is not the case. But there 
is a close relationship with a triple coalescence. There 
are two exceptional points located in the circular area of 
the parameter space circle of the type (IT4|) used in Fig. 
[7Ja). This can directly be shown if one chooses two dif- 
ferent parameter loops with two different center points 
(70, /o) and smaller radii S, which is done in Figs. El^b) 
and (c). Then one observes that there are two excep- 
tional points at which two of the three resonances form 
a square root branch point. The resonance denoted by 
down triangles in Fig. [7] is involved in both exceptional 
points. For the large parameter space circle used for Fig. 
Wia) the two different exceptional points cannot be re- 
solved and the eigenvalue paths form the permutation 
of three resonances which looks like a triple coalescence 
in the form of a cubic root branch point. Although no 
exact triple coalescence was found the finding of the sit- 
uation depicted here is already a remarkable result. Two 
exceptional points located close to each other in the pa- 
rameter space are required. Furthermore, only three res- 
onances may be connected with these exceptional points, 
i.e., one of them must be connected with both branch 
points. If three additional parameters were available, a 
coalescence of both exceptional points would be possible 
or, in other words, the additional parameters would be 
necessary to shift the three resonances such that they 
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FIG. 7: (Color online) Illustration of a structure in the vicin- 
ity of two exceptional points in which three resonances are in- 
volved. Each of the resonances is denoted by a different sym- 
bol and color. The inset of Fig. (a) shows the three parameter 
space loops (circles) and the position of the two exceptional 
points (up triangles) . In Fig. (a) the parameter space loop (in- 
dicated by dots) is chosen such that both exceptional points 
are located within the circular area. The resonance marked by 
down triangles forms a branch point singularity with two fur- 
ther states at two different parameter values, which is shown 
by two additional parameter space loops (small circles in the 
inset) in whose circular area only one of the two exceptional 
points is located (Figs, (b) and (c)). The solid lines mark the 
eigenvalues of the three-dimensional model (|28b|l - (|28d|l fitted 
to the numerical data. 



form a cubic root branch point in the five-dimensional 
parameter space. 

Of course, it is not possible to describe a behavior of 
this kind with the two-dimensional matrix models used 
so far. At least a three-dimensional model is required 
to simulate three energy eigenvalues connected with each 
other. Indeed, it can be shown that it is possible to 
reconstruct the structures shown in Fig. [7] by a three- 
dimensional matrix model. Similar to the ansatz (jI5p we 
expand the matrix elements in a power series in the two 
field strengths 7 and / around a center point, and specif- 
ically assume the matrix to be symmetric. To model the 
behavior of its eigenvalues A, we fit the coefficients of the 
characteristic polynomial 



+aX^ + bX + c = 



(28a) 



to the exact numerical results. Note that one has the 
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familiar relations 



a = -(Ai + A2 + A3) , 
b = A1A2 + A1A3 + A2A3 
c = -A1A2A3 . 



(28b) 
(28c) 
(28d) 



For the discussion in Fig. [7] a power series expansion of 
the coefficients a, 6, and c up to third order in both 
field strengths was included, which leads to 30 terms 
for all three coefficients, i.e., 10 combinations of field 
strengths are required to obtain 10 triples of eigenval- 
ues for the relations (|28bp - (|28dp . The eigenvalues in 
the three-dimensional matrix model are shown as solid 
lines in Figs. [TKa), (b), and (c) and agree very well with 
the numerically exact resonances. This shows that it 
is sufficient to only take the three resonances into ac- 
count to explain their behavior. The influence of fur- 
ther resonances can be ignored in the investigation of 
the threefold permutation. The model can even be used 
to predict the positions of the two exceptional points 
located within the parameter space loop. For the case 
shown in Fig. [D^a) the model predicts the positions of 
the two exceptional points labeled 15 and 16 in Ta- 
ble U at (71 = 6.12 X 10-3, /i = 2.53 x 10-*) and 
('-)/2 — 6.15 X lO^'^, /2 = 2.68 X 10^^), respectively, which 
well approximates the results of the full quantum treat- 
ment. 



VIII. CONNECTION WITH AVOIDED LEVEL 
CROSSINGS 

There is a close relation between avoided level cross- 
ings of the real energies of bound states and exceptional 
points. As has been demonstrated, the level repulsions 
of bound states of a Hermitian Hamiltonian which de- 
pends on one real parameter are associated with an ex- 
ceptional point if the parameter is continued into the 
complex plane [H, Q . A similar effect appears with reso- 
nances of open quantum systems. Here one can observe 
crossings or avoided crossings of either the positions or 
the widths of the resonances for lines in the parameter 
space which do not run over the exceptional point. This 
behavior has, e.g., been discussed for the resonances in 
microwave cavities [TT| . 

Avoided level crossings of the energies or widths can 
also be observed in spectra of the hydrogen atom in ex- 
ternal fields. Fig. [5] shows the real part of the complex 
energy as a function of a parameter a which defines a 
straight line of the form 



7 = 0.80516 X lO^^a 
/ = 0.32169 X lO^^a 
7// = 25.03 



(29a) 
(29b) 
(29c) 



in the (7, /)-space. The line is chosen such that it passes 
directly through the exceptional point labeled 12 in Table 
m The position of the exceptional point is marked by the 
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FIG. 8: (Color online) Real part of the complex energy as a 
function of the one-dimensional parameter q, 7 = 0.80516 x 
IQ-^a, / = 0.32169 x 10"^ a, defined in Eqs. (pSal) and (l29b)l . 
The blue circle marks an exceptional point. Avoided cross- 
ings (marked by blue squares) appear, which are related to 
exceptional points. Only resonances with |Im(i5)| < 0.0005 
are drawn to not overload the figure. The black frame marks 
the region investigated in detail in Fig. Oa) to demonstrate 
how an avoided crossing leads to an exceptional point if the 
parameters are varied. 



blue circle in the figure. In the vicinity of the parameter 
space line further exceptional points are located. Since 
the line does not pass through these critical parameter 
values, the exceptional points only manifest themselves 
by avoided crossings in the energies or the widths. The 
avoided energy (real part) crossings are shown in the fig- 
ure. They are marked by blue squares. The close connec- 
tion between avoided crossings and branch point singu- 
larities of exceptional points can be shown very directly. 
In all cases marked in Fig. [5] it is possible to vary the pa- 
rameters 7 and / such that a coalescence of the two eigen- 
values forming the avoided crossing is achieved. That is, 
all avoided crossings shown in the diagram are associated 
with exceptional points of the corresponding resonances. 
In all cases where the description of the energy range 
under consideration is possible with the model i.e., 
always if only two resonances are involved, the adjust- 
ment of the parameters 7 and / leads to a branch point. 
In other words, exceptional points are indeed found to be 
the origin of narrow avoided level crossings. 

A more detailed illustration is given in Fig. [9l The 
region marked by the black frame in Fig. [5] is magni- 
fied in Fig. [DJa). Three of the avoided level crossings 
are included and, in addition, the imaginary part of the 
resonances is shown. The encounter of two resonances 
marked by the arrow labeled 1 shows an avoided cross- 
ing in the real part and a crossing in the imaginary part. 
The arrow 2 exhibits avoided crossings in the real as well 
as in the imaginary part. The encounter marked by ar- 
row 3 shows an avoided crossing of the real part which 
is accompanied by a crossing in the imaginary part. Its 
behavior changes on a neighboring line in the parameter 
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FIG. 9: (Color online) (a) Detailed illustration of the region 
marked by the black frame in Fig. [S] Three avoided en- 
ergy crossings are present. In addition, the imaginary part is 
shown, which exhibits two crossings (arrows 1 and 3) and one 
avoided crossing (arrow 2). (b) The same region is shown for 
the field strengths varied on the line 7 = 0.65378 x 10~'^q', 
/ = 0.28213 X lO-^a' given by Eqs. ifSOa)) and (|30b)) . One 
of the three exceptional points is located on this line and the 
corresponding avoided energy crossing marked by arrow 3 has 
been transformed into a branch point singularity, where both 
the real and imaginary parts are identical. The other two en- 
counters show at least one avoided crossing in the real or in 
the imaginary part. 



space defined by 



7 = 0.65378 X IQ-^a' 
f = 0.28213 X lO^^a' 
7// = 23.17. 



(30a) 
(30b) 
(30c) 



Now the avoided energy crossing is changed into a cross- 
ing, which can be seen in Fig. [9{b). The line (|30ap . 
()30bp runs exactly through the corresponding exceptional 
point. As a consequence the avoided crossing is trans- 
formed into a branch point degeneracy at which both the 
real and imaginary parts are identical, whereas the other 
two encounters do not form a degeneracy. They belong 
to different exceptional points which are not hit by the 
line. The encounter labeled by arrow 2 forms, again, an 
avoided crossing both in the real and in the imaginary 
part. The third encounter, which is marked by arrow 
1 has inverted its behavior. Now, it forms a real part 
crossing and an avoided crossing of the imaginary part. 

The connection between avoided crossings and excep- 
tional points provides an additional possibility of detect- 
ing the branch point singularities in spectra of the hy- 
drogen atom in external fields. Exceptional points can 
be found by plotting the real part of the complex reso- 
nance energies as a function of one parameter, similar to 
a in Eqs. p9a|) and (j29b|) . If avoided crossings of the en- 
ergy are found, the identification of the exceptional point 
can be performed by the method presented in Sec. IIV AI 



IX. CONCLUSION AND OUTLOOK 

Exceptional points are a feature that can emerge in 
parameter-dependent open quantum systems with decay- 
ing unbound states. They are branch point singularities 
at which two eigenstates of a non-Hermitian Hamiltonian 
coalesce. The topic of this paper was to investigate these 
exceptional points in the spectra of atoms in external 
fields. It has been shown, that the branch point singular- 
ities can be found by the permutation of two eigenvalues 
when an exceptional point is encircled in the parameter 
space. This method works reliably. 

The effects of exceptional points on spectra of the hy- 
drogen atom in external fields have been analyzed in de- 
tail. In a close region around the branch points, it is pos- 
sible to describe the two branching eigenstates by 2 x 2- 
matrix models and to explain the structure of the loops 
the eigenvalues traverse for closed paths in the param- 
eter space. The study of dipole matrix elements in the 
vicinity of branch points has revealed remarkable proper- 
ties. While single dipole matrix elements diverge in the 
presence of exceptional points this is not the case for ob- 
servable physical quantities such as the photoionization 
cross section. Both behaviors can be explained by a sim- 
ple matrix model, which provides an excellent qualitative 
description of the properties in the local vicinity of the 
branch point singularity. 

Beyond the typical square root branch points the 
spectra of the hydrogen atom exhibit structures in 
which three resonances are connected via two exceptional 
points. One of the three resonances is involved in both 
branch points. This finding is especially important be- 
cause it is closely related to a cubic root branch point, at 
which all three resonances coalesce. The direct adjust- 
ment of such a cubic root branch point would require five 
external parameters and therefore is not possible in the 
system investigated here. Furthermore, the close relation 
of avoided level crossings and exceptional points could be 
confirmed in the resonances of the hydrogen atom in ex- 
ternal fields. 

As an outlook, it should be possible to extend the two- 
dimensional matrix model introduced in this paper to 
develop a further, and possibly very fast, method to de- 
termine the exact position of an exceptional point. Once 
its existence has been detected by the permutation of two 
eigenvalues the fit to the two-dimensional matrix model 
can be used to predict the position of the degeneracy in 
parameter space by solving Eq. l|16bp for Ai — A2 which 
is a very easy task. The improved position of the excep- 
tional point can be used to perform a smaller parameter 
space circle and to apply the procedure iteratively. 

With the procedure proposed in Ref. ^4d\ it should be 
possible to verify experimentally the effects of exceptional 
points on the spectra of atoms in external fields discussed 
in this paper. In particular, the permutation behavior to 
prove their existence, the infiuence on the photoioniza- 
tion cross section, and the close relation to avoided level 
crossings should be observable. 



12 



Acknowledgments 



graduiertenforderung of the Land Baden- Wiirttemberg. 



This work was supported by Deutsche Forschungsge- 
meinschaft. H.C. is grateful for support from the Landes- 



W. D. Heiss, Eur. Phys. J. D 7, 1 (1999). 

T. Kato, Perturbation theory for linear operators [25 
(Springer, Berlin, 1966). 

W. D. Heiss and A. L. Sannino, J. Phys. A 23, 1167 
(1990). [26 
W. D. Heiss, Phys. Rev. E 61, 929 (2000). 

M. V. Berry and D. H. J. O'Dell, J. Phys. A 31, 2093 [27 
(1998). 

P. Cejnar, S. Heinze, and M. Macck, Phys. Rev. Lett. [28 
99, 100601 (2007). 

S. Klairnan and L. S. Cederbaum, Phys. Rev. A 78, [29 
062113 (2008). 

S.-Y. Lee, J.-W. Ryu, J.-B. Shim, S.-B. Lee, S. W. Kim, [30 
and K. An, Phys. Rev. A 78, 015805 (2008). 
M. Miiller and I. Rotter, J. Phys. A 41, 244018 (2008). [31 
M. Philipp, P. von Brentano, G. Pascovici, and [32 
A. Richter, Phys. Rev. E 62, 1922 (2000). 

C. Dembowski, H.-D. Graf, H. L. Harney, A. Heine, [33 

W. D. Heiss, H. Rehfeld, and A. Richter, Phys. Rev. 

Lett. 86, 787 (2001). [34; 

C. Dembowski, B. Dietz, H.-D. Graf, H. L. Harney, 

A. Heine, W. D. Heiss, and A. Richter, Phys. Rev. Lett. [35 

90, 034101 (2003). 

C. Dembowski, B. Dietz, H.-D. Graf, H. L. Harney, [36 

A. Heine, W. D. Heiss, and A. Richter, Phys. Rev. E 

69, 056216 (2004). [37 

B. Dietz, T. Friodrich, J. Metz, M. Miski-Oglu, [38 
A. Richter, F. Schafer, and C. A. Stafford, Phys. Rev. [39^ 
E 75, 027201 (2007). 

T. Stehmann, W. D. Heiss, and F. G. Scholtz, J. Phys. [40 
A 37, 7813 (2004). 

M. V. Berry, Proc. R. Soc. A 392, 45 (1984), ISSN [41 

00804630. 

N. Moiseyev, Phys. Rep. 302, 212 (1998). [42 
W. P. Reinhardt, Ann. Rev. Phys. Chem. 33, 223 (1982). 
Y. K. Ho, Phys. Rep. 99, 1 (1983). [43 

O. Latinne, N. ,L Kylstra, M. Dorr, J. Purvis, M. Terao- 
Dunseath, C. J. Joacliain, P. G. Burke, and C. J. Noble, [44 
Phys. Rev. Lett. 74, 46 (1995). 

H. J. Korsch and S. Mossmann, J. Phys. A 36, 2139 [45 
(2003). 

E. Hernandez, A. .Jauregui, and A. Mondragon, J. Phys. [46 
A 39, 10087 (2006). 

E. M. Graofe, U. Giinther, H. J. Korsch, and A. E. [47 

Niederle, J. Phys. A 41, 255206 (2008). 

P. von Brentano and M. Philipp, Phys. Lett. B 454, 171 



(1999) . 

M. K. Oberthaler, R. Abfalterer, S. Bernet, J. Schmied- 
mayer, and A. Zeihnger, Phys. Rev. Lett. 77, 4980 
(1996). 

H. Cartarius, J. Main, and G. Wunner, Phys. Rev. A 77, 
013618 (2008). 

K. Rapedius and H. J. Korsch, J. Phys. B 42, 044005 
(2009). 

A. L. Shuvalov and N. H. Scott, Acta Mech. 140, 1 

(2000) . 

S. Klaiman, U. Giinther, and N. Moiseyev, Phys. Rev. 
Lett. 101, 080402 (2008). 

S. Pancharatnam, Proc. Indian Acad. Sci. A 42, 86 
(1955). 

M. V. Berry, Current Science 67, 220 (1994). 

J. Wiersig, S. W. Kim, and M. Hentschel, Phys. Rev. A 

78, 053809 (2008). 

S. Freund, R. Ubert, E. Flothmann, K. Welge, D. M. 
Wang, and J. B. Delos, Phys. Rev. A 65, 053408 (2002). 
J. Rao, D. Delande, and K. T. Taylor, J. Phys. B 34, 

L391 (2001). 

T. Bartsch, J. Main, and G. Wunner, Phys. Rev. A 67, 

063410 (2003). 

T. Bartsch, J. Main, and G. Wunner, Phys. Rev. A 67, 

063411 (2003). 

J. Main and G. Wunner, Phys. Rev. Lett. 69, 586 (1992). 
.J. Main and G. Wunner, J. Phys. B 27, 2835 (1994). 

G. Stania and H. Walther, Phys. Rev. Lett. 95, 194101 
(2005). 

H. Cartarius, ,J. Main, and G. Wunner, Phys. Rev. Lett. 
99, 173003 (2007). 

U. Giinther, I. Rotter, and B. F. Samsonov, J. Phys. A 
40, 8815 (2007). 

A. P. Seyranian, O. N. Kirillov, and A. A. Mailybaev, J. 
Phys. A 38, 1723 (2005). 

P. Schmelcher and L. S. Cederbaum, Chem. Phys. Lett. 
208, 548 (1993). 

D. Delande, A. Bommier, and J. C. Gay, Phys. Rev. Lett. 
66, 141 (1991). 

R. B. Lehoucq, D. C. Sorenscn, and C. Yang, ARPACK 

Users' Guide (Siam, Philadelphia, 1998). 

T. N. Rescigno and V. McKoy, Phys. Rev. A 12, 522 

(1975). 

W. D. Heiss, J. Phys. A 41, 244010 (2008). 



